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Abstract —We consider a dense urban cellular network where 
the base stations (BSs) are stacked vertically as well as extending 
infinitely in the horizontal plane, resulting in a greater than two 
dimensional (2D) deployment. Using a dual-slope path loss model 
that is well supported empirically, we extend recent 2D coverage 
probability and potential throughput results to 3 dimensions. We 
prove that the “critical close-in path loss exponent” ao where 
SINR eventually decays to zero is equal to the dimensionality d, 
i.e. ao < 3 results in an eventual SINR of 0 in a 3D network. We 
also show that the potential (i.e. best case) aggregate throughput 
decays to zero for ao < d/2. Both of these scaling results also 
hold for the more realistic case that we term SD"^, where there 
are no BSs below the user, as in a dense urban network with the 
user on or near the ground. 


Contributions: The contribution of this letter is to extend m 
to three dimensions. The typical 3 D case corresponds to a user 
sufficiently high off the ground in a dense urban environment 
that they see an appreciable number of BSs in every direction. 
We also consider a case we term 3 D+, which is a special case 
of 3 D with BSs extending overhead in the positive direction 
only, corresponding more closely to a user on the ground. 
We compute the probability of coverage (SINR distribution) 
and the potential throughput for both of these cases. Then we 
compute the critical values of the close-in path loss exponent 
for which SINR and throughput respectively go to zero as the 
density goes to infinity for both the 3 D and 3 D+ cases. 


I. Introduction 

Cellular networks have been continually densifying their 
base stations (BSs) since their inception, driving most of the 
increased throughput over the past several decades Q. The 
urban environments where spectrum is most scarce and high 
capacity most critical are themselves continually densifying, 
particularly vertically, albeit at a slower rate. Future cellular 
deployments, featuring small low power base stations, will 
extend in the vertical direction as well as in the 2 D plane. 
Recently, El showed that for a 2 D network, the coverage 
probability for a given SINR value decays to zero as the 
network is heavily densified if the close-in (i.e. within a 
distance Rc) path loss exponent ao is less than 2, regardless 
of the other path loss exponents outside Rc- The potential 
throughput, i.e. the best case aggregate throughput, still grows 
at-least sub-linearly if ao > 1. 

Other Related Work: The coverage probability of cellular 
wireless systems has been studied in detail in recent years 
using stochastic geometry 0, a for various 2 D deployment 
scenarios. In 0 , a general d dimensional Poisson Point Pro¬ 
cess (PPP) BS deployment is considered and an equivalent one 
dimensional (ID) PPP is derived to compute the SINR and rate 
coverage for highest instantaneous power based association. 
The dual-slope path loss model 0 is a generalization of 
standard single slope path loss model, which was analyzed 
in 0 and is the focus of this letter. It is well supported by 
many measurements 0, 0, 0 and shown to be very close 
to many scenarios of current interest including indoor uni, 
LTE El], E2l and millimeter wave ES, HI- 
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II. System Model 


We consider a downlink cellular network with BSs located 
in a (J-dimensional space according to a Poisson Point Process 
(PPP) $ = {xi : Xi € R'*} with intensity A. The average 
number of BSs in a d-ball of radius r located at origin for 
such a PPP is given by 

A(i3(r)) = Vy\ 


where VdT'^ is the volume of d-ball B{r) and 14 is a constant 
dependent on the dimension d. Note that V2 = tt and V3 = |7r, 
while for the 3 D+ case, BSs are located only in a half sphere, 
so ^30+ = V3/2 = I?!. With a slight abuse of notation, we 
will use c? = 3 D and 3 D+ to differentiate between 3 D and 
3 D+ cases where necessary, keeping in mind that d = 3 for 
both cases. 

We assume a dual-slope path loss model which has two 
different path loss exponents, as in 0, given as 


r ““, for r < Rc 
^ for r > Rc^ 


where Rc is the critical distance, ao is the close-in path loss 
exponent and ai is the long-range path loss exponent, with 
T] = “““ a constant to provide continuity. We require 0 < 

cto < cti, and assume Rayleigh fading for all links, therefore 
the received power at the origin from the BS located at Xi 
is given by Pi = hi£{\\xi\\), where hi’s are i.i.d. exponential 
random variables with mean 1. 

We assume that the user connects to the BS providing the 
highest average received power (i.e. closest BS) and denote 
this BS by index 0 . Therefore the SINR is 


SINR = 


_^of(lkoll)_ 

+EiG$\{o} 


where is the noise variance. 
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We are interested in the following two performance metrics. 

Definition 1, The downlink coverage probability in d dimen¬ 
sions is 

PsiNR(A,r)=P[SINR>T], 
which is equivalently the ccdf of the SINR. 

Definition 2. The potential throughput rf' captures the average 
number of bits that can be transmitted per unit area per unit 
time per unit bandwidth, assuming all BSs transmit (i.e. full 
buffer model), and is 

rf(A,T) = log2(l + r)AP^ir,R(A,T). 

It has units of area spectral efficiency: bps/Hz/m^. 


III. Probability of Coverage 
In this section, we first derive the coverage probability 
expression for general path loss function and present the 
simplified expression for the dual-slope path loss model. 


Lemma 1. The coverage probability with a general path loss 
function is 


Psinr(A,T) = AF, / 
^0 


/ 


X exp 


/ 


-XVdy 


V 


1 + 


V 


rp _|_ 


-—dt \dy. (1) 


mv)^) 


For the 3D and 3D+ cases, Vd is ^ and respectively. 


Proof: See Appendix lAl ■ 

We observe that expression for 3D and 3D+ are the same, 
except for the value of Vd- Also, the effective BS density for 
3D+ is A when compared to the 3D case. Thus, from now 
on, we will consider only the 3D case for analysis with the 
understanding that such results can be trivially converted to 
3D+. Fig. [T] validates Lemma [T] by comparing a simulation 
of the system model with the analytic expression given in 
the Lemma. It also shows that the SINR coverage of a 3D+ 
deployment with density A is equal to that of a 3D deployment 
with density A. 

Lemma [T] can be further simplified for the dual-slope case 
to give the following Theorem. 


Theorem 1. The downlink coverage probability for a general 
d-dimensional PPP BS deployment under the dual-slope model 
is given as 


’sinr(A, 7") — AVdi?^ / e 


d f -\VaRil(T,r)-T(7^ 


dr 


/ OO q: 


where 


I{T,r) =C 

-rC 


ao 1 




d'Tr^ 

22 . i. 

d ’ T 


Oil 2R 


r - 1, 


C {b, z) — 2 F 1 ( Ij Ti 1 + 



SINR Threshold (dB) 


Fig. I. SINR coverage probability for 3D+ and 3D BS deployment with 
ao = 3.3, ai = 5 , cr^ = = 0 . 4 . 


with 2 F 1 (a, b, c, z) being the Gauss hypergeometric function. 

Proof: See Appendix iBl ■ 

Before going further, we will also compute the SIR coverage 
probability assuming noise to be zero which mimics the inter¬ 
ference limited case. SIR coverage probability tightly upper 
bounds SINR coverage probability for dense deployments and 
is given as 


PsiR(A,r) 

= WdRi -f ^ 

= XVdRi / . 

Vo 


The SNR coverage probability can be found similarly by 
letting the interference go to zero. 

The following Lemma establishes the relationship between 
the 2D case considered in ||2l and the general d-dimensional 
case. 


Lemma 2. The probability o/fSIR, SINR and SNR) coverage 
for a general d dimension PPP BS deployment with parame¬ 
ters ctoi ctii A, i?c, A equal to the probability of coverage 
for a 2D system with af,, a'l, X\Rc, if 


Qn — 




a, = -ai, 




/ /\2 2 r>o‘o — ao 

[a ) = a Rc 


Proof: Proven easily by substituting the respective param¬ 
eters in (Ell for d = 2 and observing the exact same expression. 


Following the similarity of SINR expression to that in EIj 
it can be shown that El Lemma 2 ] and El Theorem 2 ] will 
also be valid for the general d dimensional case. Building on 
these results and Lemma |2l we state the following Theorem. 
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Fig. 2. SINR and SIR coverage scaling vs. network density (A) for 3D 
deployment, with ao = [2.5, 3.5], ai = 4, Rc = 0.4, tr^ = 1. 



Fig. 3. Potential throughput (t;) scaling with network density for 3D 
deployment. Here, oq = [1,1-5, 2, 3], oi = 4, Rc = 0.4, = 1, T = 1. 


Theorem 2. Under the dual-slope path-loss model, the SIR 
and SINR coverage probability of a general d-dimensional 
system go to 0 as X ^ oo for ao < d. 

Proof: See Appendix O ■ 

The above Theorem is true for general d dimensional 
deployments and hence is valid for both the 3D+ and 3D cases. 
It provides the critical values of the close-in path loss exponent 
below which the coverage probability goes to zero. Theorem 
implies that for both the 3D and 3D+ scenario, the critical 
value of ao is 3. It is very common for the path loss exponent 
of short range systems to be less than these ao values, so 
this is seemingly an important concern for future ultra dense 
networks. 

Fig. |2] shows the behavior of SINR and SIR coverage 
probability (P|inr and PHr) for a 3D BS deployment as 
the network density varies. It can be observed that for all path 
loss exponents, P|inr first increases as A increases. After a 
critical limit of A, P^^r starts decreasing. For ao less than 3, 
Pg^R goes to zero while for ao < 3, P^^r asymptotically 
becomes a nonzero constant as A —)• cxd. For lower A, Pg^ 
corresponds to coverage probability for single slope path loss 
model with a\ and for higher A, Pg^ corresponds to that with 
ao. We can observe that Pg^ goes to zero for ao < 3 as A 
goes to infinity. 

IV. Potential Throughput 

We now turn to the potential throughput scaling with 
density. 

Theorem 3. Under the dual-slope model, as X ^ oo, the 
potential throughput rf 

1 ) grows linearly with X if ao > d, 

2 ) grows sublinearly with rate if I <a<d, 

3) decays to zero if ao < |. 

Proof: Using Lemma |2] we can prove that the potential 
throughout in a general d dimensional BS deployment is 


connected to that of 2D case by the following relation: 


d(\ T-' ^ „ D ,P2D f \l rp 


^ 0 ’ 


a'i,i?c,(cr')^)- 


Using the above relation and ||2] Theorem 3], all three results 
of the Theorem |3] can be easily proven. For ao > d, E 
Theorem 3] states that potential throughput in 2D case scales 
linearly with A' for a^ > 2. Therefore, the potential throughput 
in the general d-dimensional case will scale linearly with A if 
a = a' d > d. Similarly, the other two results for ^ < ao < d 
and ao < f can be obtained. ■ 

Fig.[3]shows the scaling of potential throughput with respect 
to A for a 3D deployment. As expected, for 1.5 < ao < 
3 the potential throughput scales only sub-linearly. Theorem 
|3] provides a theoretical basis for understanding the gain in 
throughput vs. the cost of densification. 

We conclude by noting that the SINR throughput and SINR 
scaling results can be easily extended to more than two path 
loss exponents. Owing to the equivalency between the 2D case 
discussed in E and the general d dimensional case. Theorem 
12] and Theorem [3] can be shown to be true for the multi-slope 
path loss model also. 


Appendix A 
Proof of Lemma[T] 


Let us denote the sum interference at origin by I . Now the 
SINR coverage probability can be written as 


SINR 


= P [SINR > T] =: 


ho^iWxoW) 


> T 


g-T.V^(U£^(T/£(x))/||,„||(a:)dx (3) 


where /||a:o||(a:) is the probability distribution of the dis¬ 
tance Ijxoll of the closest (serving) BS from origin given 
as /||a;o||(a:) = VddXx‘^~^e~^^'‘^ and Ci{s) is the Laplace 
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transform of interference which can be derived as 


£ 7 ( 5 ) = E[e-"^] =E 




, M(||z 


= exp ( —A 


s£{z) 


Vddz’^ 1 . 


1 + s£[z) 

Substituting the values of /||a;o||(^) Ci{s) in (O, we get 


PsiNR — ^ddX / 
exp ( —X 

(a) 




£{z)Tl(.{x) 


Vddz‘^ ^dz ) dx 


1 + l{z)Tli(x) 

( \ /“^ 

VddX / 

Jo 


exp ( —XVddx° 


li i{x)/£{ux)+T 


u‘^-^du 1 dx 


where (a) is due to the substitution zjx ^ u. Now using the 
substitutions x"^ —5- y and u'^ —> f, we get O- 


Appendix B 
Proof of Theorem[T] 

Using the substitution, r = y/R^ in ([T]i, we get 


PsiNR(A,r) = Au,i?; 


7 ' 


( 


X exp —XVdR^r 


1 + 


I 


,-Ta^/i{Rar^) 


T 


Jl _|_ e(RarZ) 

\ e{R^{tr)i) 


-dt dr. (4) 


Outer integration from r = 0 to 00 in (01 can be divided 
into integration over the following two intervals: 

1) In interval [0,1); the inner integral can be written as 


1 T + t^ 

= 1 + -C 


-dt 


T 


4 T + 


-df 


Do 1 


-C 


ctp j_ 
d ’ T 


r \d' Tr^ . 

+ icf-^,-Tr^r-^)-7. (5) 

r \ a Jr 

2 ) In interval [ 1 , oo): the inner integral can be written as 

n oo rj^ 


d( = c(-^,T) 


- 1 . 


J1 T 1 'i 
Using (0 and (| 6 ]l in (01i, we get (|2]i. 


( 6 ) 


Appendix C 
Proof of Theorem[2] 

This proof is similar to the proof of Proposition 1 in 121. 
As the SINK coverage probability is always less than SIR 
coverage probability, it suffices to show the proof for SIR only. 
Using Lemma 1 and taking cr^ = 0, we can upper bound the 


SIR coverage probability as following: Psir(Ai T) 


< XVd j exp I -WdV I 1 + 


I 


max (1, 


-dt 


T + 


JM) 

my)h 


dy 


XVd 




= XVdRt J exp I -XVdRiy I 1 + 


T 


\ T + t^ 


-dt dy 


^-WdRi 


The second term goes to zero as A —0. To prove the same for 
first term, consider an increasing sequence {A„}, and define 


Jn[x) = A„exp -XnVdRcV 


/i T + t-x 


It is clear that /n(x) 

fu{x) <g{x) = 


0 point wise for each x in (0,1). Also, 
1 


V,R{,h+{-^ 


dt 


g{x) is integrable on (0,1) for 0 < ag < d. So by the 
dominance convergence theorem, first term, and hence the sum 
also goes to zero which proves the Theorem. 
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